Orbital magnetic properties of quantum dots: the role of electron-electron interactions 
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We study the magnetic orbital response of a system of N interacting electrons confined in a 
two-dimensional geometry and subjected to a perpendicular magnetic field in the finite tempera- 
ture Hartree-Fock approximation. The electron-electron interaction is modelled by a short-range 
Yukawa-type potential. We calculate the ground state energy, magnetization, and magnetic sus- 
ceptibility as a function of the temperature, the potential range, and the magnetic field. We show 
that the amplitude and period of oscillations in the magnetic susceptibility are strongly affected by 
the electron-electron interaction as evidenced in experimental results. The zero-field susceptibility 
displays both paramagnetic and diamagnetic phases as a function of temperature and the number 
of confined electrons. 

PACS numbers: 73.21. La, 75.75.+a 



I. INTRODUCTION 

Much effort is currently devoted to the study of 
electron-electron interaction effects on the ground state 
and transport properties of mesoscopic and nanostruc- 
tured devicespi* 2 ^ The theoretical analysis of small quan- 
tum dots, with N < 8 electrons, costumarily employs 
sophisticated numerical tools^ Those become computa- 
tionally prohibitive for larger dots, that call for a more 
schematic approach, such as mean- field and/or semiclas- 
sical approximations. In this paper we address the spe- 
cific question of how the electron-electron interactions 
change the magnetic susceptibility of quantum dots us- 
ing the self-consistent Hartree-Fock approximation. 

Our motivation stems mainly from the intriguing ex- 
perimental data collected in the early 90's by Levy and 
collaborators^ who measured the orbital magnetic sus- 
ceptibility of an array of mesoscopic squares lithografi- 
cally inscribed in a GaAs/AlGaAs heterostructure. The 
authors found a surprisingly large paramagnetic suscep- 
tibility and a power law dependence of the zero-field 
susceptibility Xo(B) = x(B — Q>^) with temperature. 
The first result is understood by means of a semiclas- 
sical single-particle analysis, showing that the magnetic 
susceptibility in ballistic devices is determined by the en- 
closed area and the stability of the shortest classical peri- 
odic orbits of the cavitypliL&SiiSiii More specifically, the 
cavity geometry determines the classical periodic orbits 
relevant for x(_B,T), provided B is weak. For generic 
integrable systems such orbits come in families and the 
magnitude x(B,T) scales as kph, where kp is the wave 
number at the Fermi energy and L is the typical length 
scale of the cavity. In contrast, periodic orbits of chaotic 
systems are isolated and x(i?,T) oc (kpL) 1 / 2 . The semi- 
classical theoretical analysis also predicts an exponential 
decay for Xo(T), which conflicts with the experimental 
data that display a much weaker temperature depen- 
dence. The experimental results^ raised natural and im- 
portant questions, concerning the role of disorder and 



interactions. These issues triggered interesting advances 
in the semiclassical approach. The disorder in the de- 
vices we are interested in is weak and long ranged, and 
hence dominated by small angle scattering. It was con- 
vincingly shown— that such kind of disorder has little in- 
fluence on the semiclassical results derived for perfectly 
clean systems. A more far reaching and fundamental 
question is the role electron-electron interactions. In a 
sophisticated semiclassical approximation to the many- 
body problem, Ullmo and collaborators^ showed that 
interactions make x{B,T) scale as kpL for integrable 
systems and kpL / \n(kpL) for chaotic ones. There two 
key approximations that allow for a full analytical treat- 
ment presented in Ref. Il'il One approximation is to 
consider a zero-range residual interaction, that beclouds 
the exchange interaction, and the other is to give up self 
consistency. 

Further motivation stems from a more recent ex- 
periment on the orbital magnetization of quantum-dot 
arraysi* that found a magnetization value two orders of 
magnitude larger than the one predicted by the noninter- 
acting single-particle picture. Existing theoretical predic- 
tions based on the direct numerical diagonalization of the 
Hamiltonian for a few elect ronsi^iSiiLi^ and on a mean 
field approximation of the many-body problem^ indi- 
cated that the electron-electron interaction plays domi- 
nant role in the magnetic properties of quantum dots as 
clearly observed in the experiment^ 

In this work, we resume the discussion on the influ- 
ence of electron-electron interactions by studying the or- 
bital magnetic response of a A-electron interacting sys- 
tem in a self-consistent Hartree-Fock (SCHF) approxi- 
mation. The screened Coulomb electron-electron inter- 
action is modelled by the short-range Yukawa potential 
V{r) — Voe~ Kr jr. Exact results are only known for the 
N = 2 caseiiii^ The mean field approximation allows us 
to study systems up to N ~ 40 interacting electrons. Al- 
though the typical experimental dots in Refs. 00] have 
a much larger number of electrons, our study provides a 
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qualitative understanding of the interaction effects in the 
magnetic properties of the system. 

We calculate the ground state energy E g , magnetiza- 
tion, and magnetic susceptibility as a function of the rel- 
evant parameters of the system, namely, the temperature 
T, the potential range n and the magnetic field B. Our 
results show discontinuities in E g (B) at T = 0, as pre- 
viously reported^ These features are smoothed out for 
finite temperature. 

We show that the magnetic response of a two- 
dimensional electronic cavity depends strongly on the 
electron-electron interaction and the magnetic suscepti- 
bility shows an oscillatory behavior similar to de Haas- 
von Alphen oscillations in metals. Furthermore, the am- 
plitude and period of such oscillations are modified by 
the electron-electron interaction. We discriminate the 
kinetic, direct, and exchange contributions to the total 
magnetic susceptibility. The direct and exchange con- 
tributions also oscillate but with a different phase com- 
pared to the kinetic contribution. We also find that the 
susceptibility dependence on temperature and interac- 
tion strength displays non-universal features which are 
strongly dependent on the parameters. A slight enhance- 
ment of the zero-field susceptibility is seen as the num- 
ber of electrons increases. The interaction-induced mag- 
netic susceptibility shows paramagnetic and diamagnetic 
phases both as a function of temperature and interaction 
strength. A rather unexpected result is the behavior of 
the exchange interaction contribution to x(-B), which be- 
comes larger than the direct contribution as the interac- 
tion strength increases. 

The paper is organized as follows. In Sec. [H] we present 
the model system considered in this study and the mean 
field solution in the SCHF approximation. The main re- 
sults are shown in Sec. III. Section IV brings the final 
remarks and conclusions. We also include an appendix 
where we present some specific details of the SCHF nu- 
merical implementation in the presence of discrete sym- 
metries. 



II. THE MODEL 

We consider the problem of N two-dimensional (2D) 
interacting electrons in a confining potential subjected to 
an external magnetic field B perpendicular to the elec- 
tron system. Since we study the orbital contribution to 
the magnetization, we are allowed to simplify the prob- 
lem and treat the electrons as spinless. The model Hamil- 
tonian reads as 



H 



N 



N 



h(r n ) + } u(r n ,r n /) , 



n<n' 



(i) 



where r„ indicates the position of the nth electron. The 
single-particle Hamiltonian h is given by 



Mr) 



f 



2m* 



A(r) + u(r) . 



(2) 



where m* is the electron effective mass. The vector 
potential A is chosen in the symmetric gauge, namely, 
A = (— By/2, Bx/2, 0). The magnetic field is expressed 
in units of $/$o ; where $ = BA is the magnetic flux 
through the system area A and $o = hc/e is the unit 
quantum flux. We choose the confining potential u(r) 
as the 2D square well of side L that closely models the 
experiment^ 

To account for screening effects the electron-electron 
interaction v is modelled byAS 



v(r,r') = 



-re|r- r'| 



47re e r |r-r'| 



(3) 



where k gives the effective interaction range and e r is the 
background dielectric constant. For k — there is no 
screening and the bare Coulomb interaction is recovered. 
Even though the v(r, r') in Eq. is different from the 
screened electron-electron interaction in the 2DEG when 
effects of the finite layer-thickness and image charges are 
taken into account the Yukawa potential captures the 
main features of a more realistic interaction and has the 
advantage of being computationally more amenable to 
handled We recall that the semiclassical approach is 
forced to use zero-range residual interaction in order to 
make the calculations feasible. By doing so, one looses a 
handle on the exchange interaction. 

For a square dot of side i, the potential energy 
scales with 1/L while the kinetic energy scales with 
1/L 2 . Hence, as L is increased, the potential energy 
becomes increasingly more important. It is then useful 
to introduce 2 ^ an "effective strength" parameter L/a B , 
where a* B = ft 2 (47reo£ r )/m*e 2 is the effective Bohr ra- 
dius. The standard dimensionless parameter that quanti- 
fies the ratio between electronic potential and kinetic en- 
ergies is r s , that in 2D reads r 2 — A/(Nir[a B ] 2 ). Hence, 
L/a* B and r s are related as r s = (L / a* B ) / VnN . We con- 
sider here a range of parameters such that 1.5 < r s < 2, 
within the typical values of r s in the experiments^^ 

We calculate the ground-state energy in the SCHF ap- 
proximation for finite temperatures. The SCHF equa- 
tions read a a<$M$i 



h(r)<t>i(r) +J2 n >J *V*(r>(r,r')<MO 



<Mr) 



E 



n 3 /dr>*(r')«(r,r')^(r)^(r') 



ef F Mr), (4) 



where the sums run over all HF orbitals. Here = 
{exp[(ef F — //)//cbT] + I} _1 is the Fermi occupation num- 
ber of the ith HF orbital with corresponding wave func- 
tion <f>i(r) and energy £^ F . As standard, the chemical 
potential /i is determined by requiring that N = £^ n,. 
We truncate the number of orbitals and take only the 
M w 27V lowest energy states into account. 
The SCHF ground state energy is given by 



EfF 



T 



HF 



HF 



v, 



HF 



3 



i i,j 

-(<t>i<f>j\v\4>j4>i)), (5) 



where the \<fii) are the HF orbitals, self-consistent solu- 
tions of Eq. O and T HF , Vf F and V^ F are the kinetic, 
direct and exchange contributions to the ground state 
energy respectively. Actually, it is numerically less ex- 



pensive to compute Ef F as 



= - 



-HF 



2 ^ 



(6) 



The calculations of the matrix elements were done nu- 
merically, with no further approximations other than set- 
ting the numerical precision. The number of elements 
grows as M 4 , where M is the basis size. For a typical 
calculation where we take M ~ 50, a numerical evalua- 
tion of about 10 7 four-dimensional integrals is required, 
a task far from trivial. To reduce that number, we have 
used an appropriate symmetrized basis on which the HF 
potential takes a block-diagonal form. Furthermore, we 
use properties of the the Yukawa-like e-e interaction to 
reduce the four-dimensional integrals to a series of one- 
dimensional integrals on the relative polar angle. These 
steps are described in Ref. Details on the HF numeri- 
cal implementation in the presence of discrete symmetries 
are found in the Appendix. 



III. SINGLE-PARTICLE PROPERTIES 

The HF single-particle spectral properties play a ma- 
jor role for determining the system magnetization and 
magnetic susceptibility. In this section we present a gen- 
eral discussion of the model Hamiltonian single-particle 
spectrum that serves as a guide to interpret the magnetic 
results that follow. 

Figure ^ shows the single-particle HF spectrum as a 
function of the magnetic flux Q/&o for N — 10 electrons 
with r s = 1.22 for T = (left) and T — A (right). 
Here and throughout the manuscript A is the single- 
particle mean level spacing. Energy is given in units of 
h 2 /(m*L 2 ). 

A very interesting feature displayed in Fig. ^ are the 
jumps in ef F {B) for the T = case. These discontinu- 
ities appear for other values of N as well and are related 
to sudden anti-crossings at Fermi energy, i.e., involving 
the levels e^ F and £jv+i- For T = the mean field selects 
the N lowest states. Hence, by parametrically moving 
through a level crossing, the character of the last occu- 
pied state is suddenly changed, and so is the mean-field. 
Narrow crossings, with a gap Se <C A, that lead to sudden 
and very strong changes in the last occupied mean-field 
state, are responsible for the jumps. 

The jumps disappear already for very low tempera- 
tures of the order of fc^T w Se <C A, as illustrated in 
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FIG. 1: Single-particle Hartree-Fock energy levels e t as a 
function of the magnetic flux "l?/$o for (a) T — and (b) 
k B T = A. Here N = 10 and r s = 1.22. 



Fig. Hfb) . For this reason, they are of very limited rel- 
evance for the experiments we are interested in, where 
ksT 3> A. However, Fermi energy anti-crossings ex- 
plain similar features reported in the study of ground- 
state properties of quantum dots in the Coulomb block- 
ade regime, that remained so far not understood^ Al- 
beit also blurred by temperature, the later case deals 
with transport where the ground state many-body wave 
function matters and thus the Anderson orthogonality 
catastrophe can come into playi^SiS 7 . 



IV. GROUND-STATE PROPERTIES 

The magnetization m(B) and the magnetic suscepti- 
bility x{B) of an electronic cavity are obtained from 



m(B) = - 



an 

dB 



and 



X(B) = 



i d 2 n 

J^dB^' 



(7) 



In turn, the grand canonical potential Q is directly com- 
puted from the ground state HF energy by f2 = E^ F — 
TS — Nfj,, where S is the entropy and [i the chemical 
potential, both functions of B. 

In the remaining of the paper, the magnetization is 
given in units of the effective Bohr magneton /i* B = 
eh/2m*c and the magnetic susceptibility is expressed 
in units of the Landau susceptibility, namely \xl\ = 
e 2 /(127rm*c 2 ). For simplicity, the non-interacting case 
is referred to as u L/a* B = 0". 
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FIG. 2: Top panel: zero-field susceptibility \o as a function 
of the particle number N for the non-interacting (thick solid 
line with filled circles) and r 3 — 1.5 (open circles) cases at 
ksT — A/2. Middle panel: the same for the kinetic term 
(open squares). Bottom panel: the same for the exchange 
contribution (open triangles) to the magnetic susceptibility. 



FIG. 3: Zero-field susceptibility xo (circles) as a function of 
interacting strength L/a* B for N = 20 electrons and UbT = 
A/2. Kinetic (squares), direct interaction (filled triangles), 
and exchange interaction (empty triangles) contributions to 
Xo are shown for comparison. 



A. Exchange contribution 

Following Eq. J5J), one can distinguish different contri- 
butions to the magnetic susceptibility, namely 

X = X kin + X d + X x (8) 

arising from the kinetic, direct, and exchange interaction 
terms of the HF ground state energy E^ F (B). As shown 
bellow, it is instructive to compare the later with the 
susceptibility computed for non-interacting electrons in 
a square cavity, ^ non - mt . 

The zero-field susceptibility xo(T) = x{B = 0, T) is 
obtained for different N up to N = 20 at k B T = A/2. 
The results are shown on Fig. [21 The noninteracting 
susceptibility ^ non - int (filled circles) oscillates with in- 
creasing amplitude as N is varied (or equivalently, as a 
function of kpL). This is in agreement with the pre- 
vious semiclassical calculations As we include the 
electron-electron interaction included, keeping r s w 1.5 
for all N, the interacting susceptibility (open circles) sen- 
sibly deviates from the noninteracting case. We find that 
such deviations are mainly caused by the exchange con- 
tribution to the susceptibility. While the kinetic contri- 
bution (squares) resembles the noninteracting situation, 
the exchange term (triangles) exhibits rather large fluctu- 
ations, which are the main responsible for the deviations 
in the total interacting susceptibility. 

The exchange contribution to the susceptibility in- 
creases with the interaction strength L /a* B and is already 
of the order of the kinetic contribution for r s « 1.5. This 
is illustrated in Fig. 01 where we chose a maximum of the 
noninteracting Xo° n_UIt (iV), namely, N — 20 and vary 
L/a* B . For low values of L/a B , the paramagnetic kinetic 
term x km (filled squares) dominates and stays almost 



constant while the total susceptibility decreases. This 
behavior of the total susceptibility is dictated by the ex- 
change term x x ? which is initially of the order of the di- 
rect contribution. As the interaction strength increases, 
both the absolute value of the direct and exchange terms 
increase. It is interesting to notice that x x overcomes x d 
for L/a* B w 5 [r s « 0.63) and is of the same order as the 
kinetic term for L/a* B « 10 (r s « 1.26). 

This is a somewhat unexpected result, since the direct 
contribution V^ F to the ground state energy is 3 to 4 
times larger than the exchange contribution V^ HF in this 
range of r s . It turns out, however, that the magnetic V^ HF 
is more sensitive to variations of the magnetic filed than 
V rf HF , giving rise to the larger exchange contribution. 

B. Magnetic field effects 

The electron-electron interaction also induces some 
nontrivial effects in the magnetic field dependence of both 
the magnetization and the magnetic susceptibility. The 
results for m(B) and x{B) m the Coulomb (k = 0) case 
are shown in Figs. 0] and [5] for different N and L/a* B . In 
this subsection we take k B T of the order of A, to wash 
out effects due to level crossings, while preserving the 
quantum effects due to long energy range spectral corre- 
lations. 

In the noninteracting case, the magnetization curves 
for different N values are very similar. For low fields, 
small oscillations arise and, as $ increases, a positive 
magnetization phase appears. When the interaction be- 
tween the electrons is included, the magnetization curves 
for N = 10 and N — 20 are quite different (see Fig. Q). 
In particular, for L/a* B = 10, the orbital magnetization 
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FIG. 4: Magnetization m(B) as a function of B — &/L 2 for 
L/a B = (circles), 2 (diamonds), 5 (filled triangles), and 8 
(open triangles). Here TV = 10 and feT = A. Inset: TV = 20 
for L/a* B — (circles) 3 (diamonds), 5 (filled triangles), and 
10 (open triangles) at k B T = A/2. 
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FIG. 6: Magnetic susceptibility x/Xl (thick solid line) as a 
function of $/$ for TV = 10 at k B T = A/2 and L/a* B = 5.0. 
The symbols correspond to the kinetic (circles), direct (filled 
triangles) and exchange contributions (open triangles). 
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FIG. 5: Magnetic susceptibility x/IX-tl as a function of the 
magnetic flux for TV = 5 (top) and TV = 10 (bottom) with 
L/a* B =0 (circles), 2 (diamonds), 4.7 (filled triangles), and 8 
(open triangles). The temperature is k B T = A. 



is about four times larger for TV = 20 then for TV = 10. 
This suggests that large systems can display strong or- 
bital magnetization effects, in line with the experimental 
results^ 

For TV = 5, a clear oscillatory behavior is seen on x(T3) 
(see top of Fig. |3J. The amplitudes are of order \ and, 
as the interaction strength is increased, both the ampli- 
tude and frequency also increase. This behavior is similar 
to the de Hass-von Alphen effect observed in metals. In 
that case, both the amplitude and frequency of x{B) os ~ 
dilations are proportional to the chemical potential of 



the system. In the single-particle effective potential ap- 
proximation, an increase in the interaction strength is 
equivalent to an increase in the effective chemical poten- 
tial of the HF levels, which corroborates the analogy. 

The oscillations are modified with increasing number of 
particles, see bottom panel of Fig. [5] For larger values of 
L/a* B , x(B) displays fluctuations near $/$o = 10 due to 
crossings between the Hartree-Fock single-particle states 
at the Fermi level, as discussed in Sec. IIIII Notice that 
although the level crossing jumps in the single-particle 
spectrum seem to have already disappeared at fc^T = A, 
see Fig. ^ they are greatly enhanced in x(T3), due the 
second derivative in Eq. Q. 

We observe that the exchange part plays an important 
role to for the computation of the B dependence on the 
magnetic susceptibility x- bi Fig. all contributions are 
shown as a function of $/$o for a lower temperature, 
namely ksT = A/2 and r s = 0.9. The pronounced para- 
magnetic peak around $/$o ~ 10 (also observed at the 
bottom panel of Fig. [SJ) is mainly due to the exchange 
contribution x x - Notice that frequently % d and x x give 
contributions of opposite signs appears as $/$o is varied. 

It is worth mentioning that these results are somewhat 
different from the ones obtained by the exact diagonaliza- 
tion for the two-electron case^ where the electron spin 
plays an important role on the orbital properties of the 
system. In that situation, singlet-triplet transitions give 
rise to low-temperature peaks of about 3|xl| of magni- 
tude in the orbital susceptibility but no oscillatory be- 
havior is found on x{B). Also, the magnetization m(B) 
in the exact TV = 2 case is always negative, even in the 
large field regime. 
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FIG. 7: Magnetic susceptibility x(B) for different interaction 
ranges (Coulomb [circles] and = L/10 [triangles]) for TV = 
5 (top) and N — 10 (bottom). 

In Fig. [7] we compare the susceptibilities in the case 
of Coulomb and short-range potentials. In order to get a 
qualitative comparison, we consider the strong-screening 
regime, using = L/10. 

In this short-range limit, the interaction between the 
electrons is exponentially suppressed for distances larger 
than the screening length k -1 and the overall magnetic 
response, which is averaged over the whole dot, is then 
similar to the noninteracting case. This can be clearly 
seen in Fig. [7Jvhere the curves for = L/10 are quite 
similar to those for L/a* B = depicted in Fig. [SJ show- 
ing that the shielding of the electron-electron interaction 
is quite effective and a noninteracting picture is a good 
approximation for the thermodynamic properties of the 
system. 

For such low screening lengths, the Fermi wavelength 
Xf ~ nL/y/N is larger than kT 1 by one order of mag- 
nitude. In the regime Xf 3> the 2D screened in- 
teraction can be well approximated by a <5-type contact 
interaction^ Our results show that the noninteracting 
picture captures the main features of such approxima- 
tion in the Hartree-Fock scheme. Therefore, one has to 
consider larger screening lengths in order to see resid- 
ual interaction effects in the magnetization and magnetic 
susceptibility. 



C. Temperature dependence 

The zero-field susceptibility decays with temperature, 
as expected, but the decay rate is rather large. For 
ksT/A = 3, xo is already negligible compared to the 
T = value (see Fig. In the intermediate range 

of kgT ~ A, non-universal features arise. As an exam- 
ple, the inset of Fig. |S] shows our result for a dot with 
N = 5 electrons: Xo(T) is negative (except for T ~ 0) 



FIG. 8: Temperature dependence of the zero-field susceptibil- 
ity Xo(T) for N = 10 and L/a* B =0 (diamonds), 2 (squares) 
and 5 (triangles). In the inset, the same is shown for N — 5 
with an additional curve for L/a* B =8 (circles). 

and has pronounced minima at ksT = A/2, even in the 
non-interacting case. For TV = 10, Xo(T) displays both 
positive (for fc^T < 0.4A and fc^T > A) and negative 
values (0.4A < ksT < A), indicating a diamagnetic- 
paramagnetic transition. In both cases, the interaction- 
induced susceptibility, X — X non_mt i is positive, indicating 
a paramagnetic contribution due to the electron-electron 
interaction. 

Such a diamagnetic-paramagnetic transition for in- 
creasing N is in agreement with previous theoretical re- 
sults. For N = 2, exact diagonalization results^ show 
that Xo(T) is diamagnetic for a temperature range up to 
10 A. On the other hand, semiclassical analysis^^ sug- 
gests a paramagnetic Xo{T) for N ~ 100 . Therefore, a 
diamagnetic-paramagnetic transition is expected in the 
intermediate particle number regime. 



V. CONCLUDING REMARKS 

We have considered the ground-state properties of a 
system of N interacting electrons confined in a 2D geom- 
etry and subjected to a perpendicular magnetic field in 
the finite temperature SCHF approximation. The mag- 
netic susceptibility was calculated as a function of the 
relevant parameters of the system (magnetic field, num- 
ber of electrons, temperature, and strength and range of 
the particle interaction). The ground-state energy was 
obtained for both the Coulomb interaction and the short- 
range Yukawa potential. 

Our results show that the electron-electron interaction 
introduces nontrivial effects in the magnetic properties of 
the system. The magnetic susceptibility shows de Hass- 
van Alphen-like oscillations which are enhanced as the 
interaction strength increases. The magnetization in- 
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creases when more electrons are added in the dot, which 
indicates that a strong orbital magnetization should be 
expected for larger systems. For a higher number of elec- 
trons, new features arise, including strong diamagnetic 
fluctuations as a function of the magnetic field. 

The zero-field susceptibility Xo(T) shows both para- 
magnetic and diamagnetic phases as a function of the 
temperature. We found that Xo(T) — > as T increases 
and the susceptibility induced by interaction is positive, 
yielding a paramagnetic contribution to \o irrespective 
of the value of N. However, non-universal ./V-dependent 
features appear in the intermediate temperature range of 
k B T » A. 
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APPENDIX A: HARTREE-FOCK EQUATIONS IN 
MATRIX FORM 

We solve the HF self-consistent equations @ iter- 
atively by diagonalization. The eigenfunctions of the 
square billiard of side length L, namely 



2 . (m-K \ . /rar 

tp a {x, y) = t sm \~ x ) sm y~L~y 



(Al) 



where a = (m, n), separated into the square four point 
symmetry classes, form the basis set {ip a }- The later 
are also eigenbasis of the operator R n /4 that rotates the 
coordinates by it /A. The eigenvalues namely c a — 

+1, — 1, +i and — i, label the square billiard symmetry 
classes. More details can be found in Ref. 0. 

The bottleneck for HF method is the calculation of the 
two-body electron interaction integrals, namely 



I a ^ 5 ee Jdv ydr'C°*(r)^*(r') V (r,r')^(r)^(r'). 

(A2) 

Here the advantage of separating the basis into symme- 
try classes comes into play. Many of the matrix elements 
la-yps are zero, depending on the single-particle symme- 
try classes involved. The "selection rule" can be summa- 
rized as follows: given the single-particle symmetry class 
of the states ip^ , ipZ , , and ipg 6 , the matrix element 
Ia-y/38 will be only nonzero i£i& 



Therefore, only states with the same two-particle sym- 
metry class are coupled by the interaction potential. 

Let us now express the HF equations in terms of the 
two-body residual interaction integrals /. The HF orbital 
wave functions read 



M 

«(r) = ^C fa ^(r). 



(A4) 



We typically truncate the basis set taking at least the 
M = 50 lowest square billiard energy states. The roman 
labels refer to the HF orbitals </>, whereas the greek ones 
to the basis set ip. The resulting matrix form of Eq. (0J 
is 



(A5) 










V af3 


= <p~ 


■ v*p, that are given by 






/dr'p(r' > r')^*(r)«(r J r')^(r) 


<9 




^r'p(r,r')C Q *(r)«(r,r')^(r') 



with the density matrix p given by 

M 



p(r,r')=^n i ^(r)0|(r'), 



(A7) 



where rii is the Fermi occupation number of the ith HF 
orbital. Notice that since p does not distinguish different 
symmetry classes, the HF mean field can effectively mix 
them. Hence, the HF potential reads 



M M 



v al = v tf3 ~ V af3 = ^2^2 D ~< 5 ( /q 7/3<5 ~ lory*?) 7 ( A8 ) 



7=1 5=1 

where, by introducing 

M 



5 ■ 



(A9) 



i=l 



c a 55 cp = Oy 55 eg . 



(A3) 



we eliminate one sum over the single-particle orbitals. 
The computation of (| A8|) is speed up by exploring the 
sparce nature of I a -ysp- 

Notice that the HF potential matrix elements will 
be, in the general case, complex numbers since the ba- 
sis itself is complex. The remaining of the numerical 
implementation is very standard. The convergence of 
the ground-state energies is obtained at iteration n if 
lE^-E^l/EKr <10-5. 
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